NASA Technical Paper 1104 


NASA 

TP 

IIOI1 

C.1 



Prediction of Ground Effects 
on Aircraft Noise 


in 

f-OAN 

afwl techni< 

KIUTLAND AFtJ, N. M 


S. Paul Pao, Alan R. Wenzel, 
and Paul B. Oncley 


JANUARY 1978 






TECH LIBRARY KAFB, NM 


TECH LIBRARY KAFB, NM 



ai3M57fl 


NASA Technical Paper 1104 


Prediction of Ground Effects 
on Aircraft Noise 


S. Paul Pao 

Langley Research Center 
Hampton, Virginia 

Alan R. Wenzel 

Institute for Computer Applications 
in Science and Engineering 
Langley Research Center 
Hampton, Virginia 

Paul B. Oncley 

MAN - Acoustics and Noise, Inc. 
Seattle, Washington 

iw\s/\ 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Office 


1978 




CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

THEORY OF PROPAGATION OF SOUND OVER GROUND 3 

EXPERIMENTS ON PROPAGATION OF SOUND OVER GROUND 5 

GROUND IMPEDANCE 8 

Methods of Measuring Impedance 8 

Impedance tube 8 

Free field 8 

Inclined track 9 

Direct pressure-velocity measurements 9 

Pressure-gradient measurements 10 

Calculation from flight noise data 10 

Ground-Impedance Data 10 

STANDARD PRACTICE IN CORRECTING FOR GROUND EFFECTS 12 

Modification of Physical Environment 12 

Flush-mounted microphones 12 

Raised microphones 13 

Gravel and similar partly controlled impedance surfaces 13 

Numerical Ground -Effects Correction 13 

Data correction with known ground impedance 13 

Miscellaneous techniques 14 

Analytic correction methods 14 

Cepstral technique 15 

RECOMMENDED PREDICTION METHOD 15 

PROBLEMS IN APPLICATION 21 

Finite Source 21 

Filter Bandwidth 21 

Partial Coherence 22 

Wind Gradients 24 

CONCLUDING REMARKS 24 

APPENDIX - RELATIVE IMPORTANCE OF PLANE WAVE APPROXIMATION 26 

REFERENCES 28 

FIGURES 32 


iii 



SUMMARY 


A unified method is recommended for predicting ground effects on aircraft 
noise. This method may be used in flyover noise predictions and in correcting 
static test-stand data to free-field conditions. The recommendation is based 
on a review of recent progress in the theory of ground effects and of the 
experimental evidence which supports this theory. This review shows that a 
surface wave, a recently discovered effect, must be included sometimes in the 
prediction method. Prediction equations are collected conveniently in a single 
section of the paper. 

Methods of measuring ground impedance and the resulting ground-impedance 
data are also reviewed because the recommended method is based on a locally 
reactive impedance boundary model. Available data support a simple model of 
ground impedance, but these data have significant scatter and it is found that 
there is need for further data, classified by terrain types such as grassland 
and desert. Experiments wherein ground impedance and the ground effects on 
noise propagation are measured simultaneously are also needed. 

Current practice of estimating ground effects are reviewed and considera- 
tion is given to practical problems in applying the recommended method. These 
problems include finite frequency -band filters, finite source dimension, wind 
and temperature gradients, and signal incoherence. 


INTRODUCTION 

The effects of propagation on noise radiated from aircraft were reviewed 
by Putnam (ref. 1) so as to develop a standard prediction method for use by the 
National Aeronautics and Space Administration. Putnam's review was one of sev- 
eral, referred to as "Key Technology Documents," which addressed different 
aspects of aircraft noise prediction. In the part of his review concerned with 
the ground effects on aircraft noise, Putnam gave concise descriptions of the 
ray acoustic analysis of sound reflections by surfaces, of sound attenuation 
over ground according to the Rudnick theory, and of the practical problems which 
are encountered in outdoor acoustic measurements. Putnam recommended empirical 
methods to predict ground effects on sound propagated from one ground station to 
another and from an aircraft to the ground. 

The present paper updates the NASA ground -effects prediction method devel- 
oped by Putnam and is designed to (1) improve the prediction methods for ground 
effects given in reference 1 , (2) include methods for correcting aircraft 
engine test-stand noise measurements to free-field conditions, (3) present a 
unified analytical prediction method for both short-distance and long-distance 
propagation over ground, and (4) define the range of elevation angle within 
which ground effects may have a significant effect on the measurement of noise 
from aircraft in flight. The prediction scheme recommended in this paper is 
intended to supplant previous empirical techniques with a unified approach. 



Significant advances have been made in recent years in understanding the 
effects of the ground surface on outdoor sound propagation. In particular, 
the rapid progress in the theory, which has occurred since 1974, has greatly 
enhanced the ability to make predictions. These advances in theory and 
experiment are summarized herein and form a basis for the recommended predic- 
tion procedures. Two recent review articles by Embleton et al. (ref. 2) and 
Piercy et al. (ref. 3) may be referred to for subjects related to ground 
effects which are beyond the scope of this paper. 

A number of situations exist in which ground effects are important in air- 
craft noise prediction and measurement. These include the prediction of com- 
munity noise resulting from aircraft operations, the correction of acoustic 
data in ground-based outdoor jet-engine noise measurement to free-field condi- 
tions, and the interpretation of data for aircraft sideline and flyover noise. 
Empirical methods have been devised to deal with each of these situations. 

The ground surface is usually modeled as an infinite flat plane on which a 
normal impedance boundary condition is prescribed. This model appears to be 
adequate for most situations of practical interest. Real media such as sand, 
soil (with or without grass cover), and snow are porous with high internal flow 
resistance and have poor wave propagation characteristics. For such media, the 
assumption of a locally reacting surface is reasonable. The strongest support 
for the validity of the impedance boundary theory comes from the cumulative 
experimental observations of ground effects in the past three decades. Indeed, 
Piercy et al. (ref. 3) have indicated that there appears to be little need for 
a more elaborate model for the description of ground effects. Therefore, most 
of the subsequent discussion herein is based on this theoretical model. 

In some situations the impedance model is inadequate. An example is 
ground covered with a thick, dense layer of vegetation (ref. 4). In this case, 
the vegetation layer should be regarded as an acoustic medium having a propaga- 
tion constant somewhat different from that of air. A strong thermal gradient 
in the air above the ground can also create an apparent condition of sound 
propagation through a layered medium (ref. 5). A layered structure can also 
occur in the solid material of the ground. In the extreme case of clay or 
hard-packed soil, the ground should be regarded as an elastic medium. A 
detailed discussion of sound propagation in layered media can be found in 
Brekhovskikh (ref. 6), and many other discussions exist in the literature for 
these special conditions. 

Perhaps the most important qualitative feature of the wave field, which 
has been brought out in the recent theoretical studies, is the surface wave. 

In contrast to the more familiar wave propagation phenomena in three- 
dimensional space (which is characterized by a more or less uniform radiation 
of energy in all directions, according to spherical spreading), the surface 
wave is effectively confined to a region near the boundary with a resultant 
decrease in amplitude with distance determined essentially by cylindrical 
spreading. An additional attenuation of the surface wave results from dissipa- 
tion of energy at the boundary; however, this additional attenuation depends on 
the value of the normal impedance of the surface. If this additional attenua- 
tion is so small that the decay of the surface-wave amplitude is determined 
primarily by the cylindrical spreading, the surface wave becomes the dominant 
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mode of propagation. Such cases are usually associated with low-frequency 
waves for real ground surfaces. 

The theories of propagation of sound over ground which have appeared since 
1974 have now effectively superseded the earlier work of Ingard, Lawhead and 
Rudnick, and others. It is desirable to suggest at this point a new prediction 
method based on the recent theoretical work. In this paper, the theories are 
reviewed first. Then a review of the existing experimental evidence of ground 
effects is given. The techniques of ground-inpedance measurement and existing 
ground-impedance data are reviewed also. 

Many currently accepted practices for ground -effects corrections are 
empirical. Although the present paper recommends a prediction method based on 
theoretical analysis, some practical limitations typical of outdoor experiments 
remain. Therefore, a brief review is given of the general status of standard 
practices in the industry. This review serves as background information for 
the application of the recommended theoretical prediction method. 

Since the present paper concentrates on the description of pure ground 
effects on sound propagation, the atmosphere above the ground is assumed to be 
homogeneous and motionless. Perturbations, such as wind, ten^jerature gradients, 
and turbulence, have not been Included in the discussions mentioned previously. 
However, variation from the ideal analytical model is certain to occur in prac- 
tical applications. Therefore, some of the common problems in application are 
discussed. 


THEORY OF PROPAGATION OF SOUND OVER GROUND 

As stated previously, the theory of sound propagation over ground is 
regarded herein as equivalent to the theory of scalar wave propagation in a 
uniform half -space, which is bounded by an infinite plane on which a constant- 
impedance boundary condition is prescribed. (The geometry of this mathematical 
model is given in fig. 1. The coordinate system, source and receiver positions, 
the direct and the reflected paths of sound propagation, and the angle of 
incidence are defined in this figure.) Only the case of radiation by a time- 
harmonic point source is considered here. Nevertheless, even within this 
rather restricted context, this review is not intended to be exhaustive. 

Instead, it is limited primarily to those contributions to the theory which the 
authors feel are most significant in terms of fundamental understanding of the 
phenomenon or else are particulau?ly useful for purposes of calculation. 

Although published theoretical work on the problem of reflection of waves 
by an impedance boundary dates back at least to the 1944 paper by Morse and 
Bolt (ref. 7), the first analytical solutions of this problem to be given in a 
form suitable for numerical calculations were those of Ingard (ref. 8) and 
Lawhead and Rudnick (ref. 9). For more than two decades the papers by Ingard 
and Lawhead and Rudnick were regarded as the standard reference works in this 
field. Recently, however, a number of in?>ortant contributions to the theory 
have appeared vAiich largely supersede the earlier work of Ingard and Lawhead 
and Rudnick. 
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Of the recent contributions, the most noteworthy are those of Wenzel 
(ref. 10), Chien and Soroka (ref. 11), and Thomasson (ref. 12). Wenzel 
obtained asymptotic results, valid when both the source and receiver are near 
the boundary, for several limiting cases determined by the values of the wave 
number, propagation distance, and boundary admittance. He also found that, 
under certain conditions, a surface wave may be a significant component of the 
total wave field. Wenzel was apparently the first to point out the existence 
and implications of the surface wave in the general case, although Brekhovskikh 
(ref. 13) had noted previously that a surface wave occurs in the special case 
in which the boundary impedance is purely imaginary."' Wenzel also noted that, 
in a certain asymptotic limit, his solution differs from Ingard's by just the 
surface-wave term. However, because these two authors used completely differ- 
ent methods, the reason for this discrepancy was not apparent. 

Chien and Soroka also obtained, as did Wenzel, asymptotic solutions for 
various limiting cases. Their results, however, are valid more generally than 
are Wenzel’s, particularly with regard to the case in which the source and 
receiver are not necessarily located near the boundary. Indeed, it appears 
that the various approximate solutions obtained by Chien and Soroka cover vir- 
tually all practical situations of interest , provided only that the field point 
is at least several wavelengths from the image source (a condition which is 
almost always met in practice). For this reason, their results are used herein 
as the basis for the recommended prediction method. 

Chien and Soroka also found that, in a certain limiting case, their solu- 
tion differs from Ingard’s by a surface-wave term, thus confirming the discrep- 
ancy noted previously by Wenzel. However, the explanation offered by Chien and 
Soroka for this discrepancy, which they ascribed to the limited range of valid- 
ity of a certain asymptotic expansion used by Ingard in the evaluation of his 
approximate solution (see ref. 11, sec. 5.2.2), reveals that these authors did 
not fully grasp the significance of the missing surface-wave term in Ingard’s 
solution, which is now known to be the result of a fundamental error in 
Ingard’s analysis. This error was brought out in a subsequent paper by 
Thomasson (ref. 12), who showed that Ingard, in the process of deforming the 
path of integration of a certain contour integral occurring in his analysis, 
failed to take proper account of a pole in the integrand. It is this pole 
which gives rise to the surface-wave term. 

The effect of Ingard’s error is to restrict the range of validity of his 
solution to essentially those values of the boundary admittance for which the 
surface-wave term does not appear. A careful examination of Rudnick’s solution 
shows that, although it was derived by a different approach than was Ingard’s, 
it is subject to a similar restriction. In contrast, the more recent results, 
since they take explicit account of the surface-wave term, are not subject to 
this restriction and, indeed, are valid for all values of the boundary 
admittance . 

'The situation in acoustics, as regards the belated recognition of the 
importance of the surface wave, thus contrasts markedly with that in electro- 
magnetic theory, in which context the surface wave has been the subject of 
much discussion, as well as considerable controversy, for over half a century. 
(See, e.g., refs. 14 and 15.) 
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In the same paper (ref. 12), Thomasson also gave an exact solutiofa (essen- 
tially a corrected version of Ingard*s exact solution) in the form of an inte- 
gral which appears to be suitable for numerical calculation. This result is 
suggested as a supplement to the approximate results of Chien and Soroka in the 
recommended prediction method. 

Thomasson (ref. 16) has also given an approximate solution, obtained from 
his exact solution, which is essentially a corrected version of Ingard's 
approximate solution. Recent contributions to the theory have also been given 
by Donato (refs. 17 and 18), Briquet and Filippi (ref. 19), and Van Moorhem 
(ref. 20). 


EXPERIMENTS ON PROPAGATION OF SOUND OVER GROUND 

Experiments regarding ground attenuation commenced as soon as the differ- 
ence between the nature of wave systems for sound propagation over a flat 
boundary and those predicted by ray acoustics or plane wave theory was recog- 
nized. Rudnick and Oncley made outdoor measurements of ground attenuation 
effects at Duke University as early as 1945. The first significant experimen- 
tal results on sound propagation over a finite-impedance boundary were those of 
Rudnick (ref. 21). He used commercial acoustic absorbing materials such as 
fiber glass in place of the ground as boundary surfaces. The experiments were 
later continued and extended by Lawhead and Rudnick (ref. 9) . In both refer- 
ences, detailed analysis accompanied the experimental observations. Among the 
more important conclusions of these studies were the following: 

(1) The sound pressure amplitude is observed to attenuate according to 
r“2 at large values of kr away from the sound source where r is the radial 
distance and k is the wave number. This observation is distinctly different 
from the expected free-field attenuation of sound pressure amplitude with dis- 
tance at the rate of r“^, according to spherical spreading. 

(2) The distance at which the rate of attenuation r“2 is attained 
depends on frequency as well as the impedance of the ground. It takes place at 
a shorter distance for a surface with low impedance than for a surface with 
high impedance. 

(3) At a given horizontal distance away from the sound source, the minimum 
received sound pressure often occurs at a finite distance above the boundary. 

Later studies by Ingard (ref. 22), Oleson and Ingard (ref. 23), and Wiener 
and Keast (ref. 24) confirmed much of these same basic phenomena of ground 
effects on sound propagation. In particular, Ingard made an inqjortant observa- 
tion with respect to the geometrical configuration tdiere both the source and 
the receiver are above the ground. He pointed out that the excess attenuation 
resulting from destructive interference between the direct and reflected sound 
paths was determined mainly by the impedance of the ground surface. The phase 
shift of ground reflection at near-grazing incidence can easily exceed 160® if 
computed according to typical values of acoustic impedance of natural ground 
surfaces. Maximum destructive interference caui be reached by an additional 
geometrical path difference of 0.05 wavelength such that the total mismatch 
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between the direct and the reflected sound waves is half a wavelength at the 
point of measurement. The ground effect in such a geometry is heavily biased 
toward destructive interference. Consequently, the peak region of the excess 
attenuation spectrum is much broader and occurs at other frequencies than those 
predicted by using ray acoustics and a boundary of infinite impedance. In ref- 
erences 22 and 24 the observed peak of excess attenuation lies between 300 and 
600 Hz. 

As is typical of many of these field experiments, the ground effect is 
often observed together with prevalent atmospheric effects on sound propaga- 
tion. The atmospheric nonuniformities usually include turbulence and often 
wind and temperature gradients. The experiments by Ingard have indicated that 
the ground attenuation effect can be studied and observed separately from most 
of the atmospheric effects. 

Some important reference data on open-field ground effects were obtained 
by Parkin and Scholes (refs. 25 and 26). Measurements were made to study the 
combined effects of the ground and wind vector on sound propagation in a hori- 
zontal direction. There were actually two groups of data. Each of these was 
taken in a different airfield where an unobstructed open field was available. 

In both cases, a jet engine of 334-kN static thrust was used as the sound 
source. Acoustic measurements were made simultaneously at eight microphones 
positioned between 20 and 1100 m away from the sound source. Approximately 
60 sets of measurements were taken in each group throughout the year. 

The data obtained by Parkin and Scholes are important for the following 
reasons: 

(1) Attenuation of sound owing to ground effect is measured for a wide 
range of horizontal distances and provides a basis for the verification of 
theoretical calculations. In particular, the data sets with zero-vector wind 
conditions are representations of the ground effects, disturbed only by turbu- 
lence but not by other meteorological effects. 

(2) Typically, data scattering occurs in field measurements. Since mea- 
surements are repeated many times over the same ground environment, the statis- 
tical confidence level of the data set as a whole is greatly enhanced. 

(3) In cases vrtiere upwind conditions are prevalent, ray acoustic theories 
predict the formation of shadow zones. However, discussions in references 24 
and 25 indicate that the observed boundary or shadow zone is frequency depen- 
dent. At low frequencies, such shadow-zone boundaries may not be observable. 
Recent theories indicate that such data may provide additional clues to the 
behavior of the surface-wave term (ref. 27). 

(4) Measurements were made at two different sites with similar arrange- 
ments and under similar weather conditions. Values of excess ground attenua- 
tion are different, vdiereas the overall characteristics remain similar. Fur- 
thermore, values of excess ground attenuation show noticeable change with 
season at each site. Hence, this data set can serve as an indirect reference 
for the characterization of ground- impedance values of grassland. 
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The determination of the acoustic impedance of typical ground surfaces is 
an important factor of the ground-effects prediction procedure. The studies of 
Delany and Bazley (refs. 28 and 29) marked the beginning of serious considera- 
tions of applying the theory to practical predictions. Based on their assump- 
tion that the acoustic impedance of soil and grassland can be modeled as typi- 
cal porous material, calculations were made to predict values of excess ground 
attenuation. The results compared favorably with the data obtained by Parkin 
and Scholes. In most comparisons the geometrical configuration dealt with 
relatively short distances with the source and receivers above the ground. 

The surface wave was not considered since the calculations were based on the 
earlier theories developed by Ingard and Rudnick. 

In references 30 and 31 a comprehensive approach was taken in the experi- 
mental observation of ground effects. Experiments were designed to observe 
both short-range (less than 20 m) and long-range (over 300 m) ground effects 
on sound propagation. The ground impedance was measured together with the 
acoustic data. Hence, theoretical predictions can be quantitatively correlated 
to the acoustic measurements. In the short-range experiments, good agreement 
is obtained between theory and experiments, regarding the amplitude and loca- 
tion of the interference patterns. However, local discrepancies at a specific 
frequency can be several decibels. In the long-distance experiment, the mea- 
surements have shown that sound attenuation owing to ground effects can be pre- 
dicted with estimated average impedance values. 

In reference 27, Piercy, Donato, and Embleton used estimated values of 
ground impedance to predict ground effects and compared them with the measure- 
ments obtained by Parkin and Scholes (refs. 25 and 26). Piercy et al. referred 
to the various components of the wave field by name. The direct wave, the 
reflected wave, and the surface wave are defined the same way as in the present 
paper. The ground wave, however, refers to a term in the asymptotic solution 
proportional to r“2. The computed results agree very well with the data of 
Parkin and Scholes. The comparison shows that the direct-wave, the reflected- 
wave, and the ground-wave components are of equal importance for sound propaga- 
tion over short distances. At longer distances, the high-frequency sound 
transmission is dominated by the combination of the direct- and reflected-wave 
components, and the received sound in low frequencies is dominated by the 
surface-wave component. Reference 30 provides an important verification for 
the correctness of the advanced theoretical work, and it stresses the impor- 
tance of the surface-wave component for sound transmission in the lower 
frequencies. 

The dominance of the direct- and reflected-wave combination at long dis- 
tances leads to an interesting conclusion. The excess ground attenuation in 
this region is a simple function of the magnitude of the complex-valued acous- 
tic impedance and the elevation angle of the sound source. In references 27, 
30, and 32 this relation was considered as a way to estimate the ground imped- 
ance by means of the excess ground attenuation measurements. However, an 
explicit explanation was not given for the relative decline of the ground-wave 
term at long distances. An analysis of the relative importance of the plane 
wave approximation is given in the appendix. 
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The previously mentioned experimental investigations represent a sequence 
of natural development in which the ground effects are observed at increasing 
levels of complexity and precision. Furthermore, the comparison between theory 
and experiment has established that the theories can accurately describe the 
physical properties of ground effects on sound propagation. Therefore, 
computer-based prediction procedure for ground effects can be established on a 
solid foundation. 


GROUND IMPEDANCE 
Methods of Measuring Impedance 

One key question relating to the study of sound propagation over ground 
which has not yet been satisfactorily answered is that of the best method of 
measuring the ground impedance. Although a number of different approaches to 
measuring the impedance of the ground have been investigated, none has emerged 
as being clearly superior under all, or even most, of the commonly occurring 
measurement conditions. 

Impedance tube .- Of the various methods which have been tried, the most 
straightforward is the ordinary impedance-tube method. In this approach, the 
impedance tube is placed vertically on the ground (sometimes the tube is driven 
several inches into the ground to provide a better seal between the end of the 
tube and the ground), and the measurements are made in the usual way. This 
technique has been applied to a sand surface by Dickinson and Doak (ref. 32) 
and to grassland by Embleton et al. (ref. 2). 

The main advantages of the impedance-tube method are the availability of 
this instrument and the simplicity of its operation. On the other hand, 
although the impedance tube would appear to be suitable for terrain such as 
grass, sand, and smooth soil, its usefulness in the case of surfaces such as 
gravel or rough soil, as well as any surface covered by a dense layer of vege- 
tation, is doubtful owing to possible difficulties in providing a proper seal 
between the surface and the end of the tube. An additional drawback, which is 
characteristic of interference techniques in general, is that the impedance- 
tube method requires accurate measurements to be made of the distance from 
the ground to the first interference minimum of the standing-wave pattern in 
the tube, with the required accuracy increasing with frequency. This require- 
ment places an upper bound on the frequency range for which this method is 
practicable. 

Free field .- In order to avoid some of the difficulties involved in apply- 
ing the impedance- tube method to real ground surfaces, Dickinson and Doak 
(ref. 32) have proposed an alternate approach, which they call the free- field 
method. The basic experimental setup consists of a sound source positioned 
above the ground together with a microphone which is allowed to move along a 
vertical axis between the source and the ground. By probing the interference 
pattern existing between the source and the ground, the ground impedance can 
be deduced. A more detailed description of the apparatus is given in 
reference 32. 
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The main advantage of the free-field method, compared to the impedance- 
tube technique, is that no tube or other wave-guiding apparatus is required, 
and hence the problem of providing a proper seal between the tube amd the 
ground surface is avoided. Since the free-field method is, however, an inter- 
ference technique, it is subject to the same frequency limitation as is the 
impedance-tube method. An additional difficulty, which plagues interference 
techniques generally, may arise in the case of surfaces such as thick grass, 
rough soil, etc., in which there is no obvious well-defined ground surface to 
use as a reference for measurements. In such cases, it is necessary to define 
a ground-surface reference level to be used as a basis for impedance measure- 
ments. All subsequent acoustic measurements, such as source and receiver 
heights, should then be based on this reference level. 

Inclined track .- A variant of the free-field method, known as the 
inclined- track method, has been used with some success by Embleton et al. 

(See ref. 2.) With this technique, impedance measurements can be made at any 
angle of incidence, instead of only at vertical incidence as with the free- 
field method. The inclined-track method thus provides a means of testing the 
hypothesis that the ground impedance is independent of incidence angle. 

Since the inclined-track method is also an interference technique, it is 
subject to the same limitations mentioned previously in connection with the 
impedance-tube and free-field methods, namely, those arising from the require- 
ments for a very accurately defined geometry and very accurate measurements, 
though these restrictions become progressively less severe in changing the 
angle from normal to grazing incidence. In order to avoid these difficulties, 
a number of alternate techniques for measuring ground impedance have been pro- 
posed. Among these are what might be called curve-fitting techniques (ref. 16, 
sec. VII; ref. 33, P- 115; and ref. 3^), the basic idea of which is to deduce 
the ground impedance by fitting calculated curves, obtained by inserting 
assumed values of the impedance into the appropriate propagation theory, to 
the obtained measurements of the sound field. Since this process involves, of 
necessity, a certain amount of trial and error, these methods are, from a 
strictly computational point of view, not as satisfying as the more systematic 
interference techniques. An additional drawback of this type of approach is 
that the experimenter frequently needs an initial estimate of the ground imped- 
ance to be sure that extraneous effects, such as the surface wave, are not con- 
taminating the measurements. 

Direct p ressure-velocity measurements .- Since specific acoustic impedance 
is, by definition, the complex ratio of sound pressure to particle velocity, 
a direct measurement of these two parameters would seem to be an obvious approach 
to impedance measurement. However, the only small microphone with a response 
proportional to velocity is a hot-wire device which is fragile and highly non- 
linear. Circuitry which linearizes the response also alters the phase so that 
measurement by this method becomes very complicated. This method has been used 
in the laboratory, usually to measure the impedance of apertures at very high 
sound pressures (see ref. 35), but it does not appear promising for outdoor 
measurements . 
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Pressure-gradien t measurements.- Although interferometric measurements are 
normally made at nodal and antinodal points for convenience, the surface imped- 
ance can actually be determined from amplitude and phase measurements at other 
points. Mechel (ref. 36) has outlined a method which appears to have some 
advantages for on-site measurements . By comparing the sound pressures and 
their phases at two small microphones, which are closely spaced on an axis nor- 
mal to the surface, Mechel has shown that the real component of the acoustic 
admittance is proportional to the gradient of phase change with height and that 
the imaginary component is proportional to the gradient of the sound pressure 
level. The gradient is directly proportional to the surface admittance if the 
measurements are made very nearly at the surface, but the admittance can be 
calculated from gradient measurements which are made some distance away. If 
the microphones are too close to the surface, their reflections alter the mea- 
sured value of surface impedance, so that this method will become more useful 
as very small microphones become available' and the accuracy of phase measure- 
ments is Improved. If the gradient is not measured very near the surface, the 
best accuracy is obtained from measurements which are made near a nodal point 
where the phase gradient is large; hence different measuring positions may be 
needed for different frequencies. 

Calculation from flight noise data.- In aircraft flyover noise measure- 
ments , acoustic data are often recorded simultaneously using a microphone which 
is mounted above the ground and a microphone which is flush-mounted on the 
ground plane. Ground -impedance information can be extracted from the ground 
interference pattern in the acoustic spectrum received by the raised microphone. 
At a given interference minimum or maximum, the magnitude of the reflection 
coefficient at the ground surface can be determined by the ratio of the mea- 
sured sound pressure amplitude to the expected sound pressure amplitude at 
free-field conditions. The frequency of the interference minimum or maximum is 
determined by the path difference and the phase shift factor of the reflection 
coefficient. Since the difference in geometrical path is known, the phase fac- 
tor of the reflection coefficient can be computed. The ground-surface imped- 
ance can be computed from the reflection coefficient provided that the angle of 
incidence is not near the zone of grazing incidence such that ray acoustic 
approximations are valid. In order to get meaningful impedance data, accurate 
position data of the source must be available and the noise data must be ana- 
lyzed on a narrow-band basis since one-third-octave spectra partly conceal the 
minimums and maximums. This method cannot reveal the entire function of imped- 
ance as a function of frequency. The effective value of ground inpedance cal- 
culated from a particular minimum may be presumed to be valid only near the 
frequency where it is measured. 


Ground -Impedance Data 

Measurements of acoustic impedance of the ground can be traced back to the 
early work by Nyborg et al. (ref. 37). However, the first systematic investi- 
gation of ground-impedance measurements appears to be the work by Dickinson and 
Doak (ref. 32). Many types of ground surfaces were examined in their work. 
These included sand, natural and tilled soil, short grass, chip granite, and 
some others. Some of the typical data are given in figure 2. One of the most 
interesting observations was the effect of moisture on the acoustic in?>edance 
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of soil. Perfectly dry sand gave the highest impedance, whereas the lowest 
impedance came when 6 to 10 percent of moisture by weight was added. Larger 
amounts of moisture increased the impedance. Dickinson speculated that small 
amounts of moisture caused grains in the soil to stick together, thereby creat- 
ing larger pores , but with larger amounts of water the pores were flooded . 

A substantial amount of ground-impedance data is given by Embleton, Piercy, 
and Olson (ref. 2). The significance of their data is that several methods have 
been used in the measurements. Cross-referencing among different sets of data 
can provide a better understanding of the variability of ground -impedance mea- 
surements. Furthermore, this data group has established a trend for the depen- 
dence of impedance as a function of frequency over an extended range. Other 
data of ground impedance are given by Thomasson (ref. 16) and Lanter (ref. 34) . 
Some data on the acoustic properties of snow are also available from Tillotson 
(ref. 38) . 

A set of data for a crushed stone field filled to a depth of approximately 
38 cm over a well-drained subsurface is given by R. H. Urban in unpublished 
correspondence to the A-21 Committee for Aircraft Noise Measurement of the 
Society of Automotive Engineers. Only the magnitude of the acoustic admittance 
is given. This set of data has been converted to the magnitude of impedance as 
a function of frequency and shown in figure 3. In the same figure, the data 
for chipped granite obtained by Dickinson et al. (ref. 32) and the data for 
grassland given by Embleton et al. (ref. 2) are given for comparison. 

Another source of ground -impedance data comes from the indirect experi- 
ments of Delany and Bazley (ref. 28). In this work, the acoustic impedance of 
porous media is obtained as a function of frequency and flow resistance of the 
porous acoustic material. By assuming that the soil is a typical porous medium, 
the acoustic impedance of soil as a function of frequency can be estimated. 
Delany and Bazley (ref. 28) used this technique to compare theoretical calcula- 
tions of ground effects with measured acoustic data, and fair agreement was 
obtained. Some comparison of the estimated ground impedance using this method 
with the data reported by Embleton, Piercy, and Olson is given in figure 4. In 
a recent study, Chessell (ref. 39) computed values of ground impedance by using 
the same model as Delainy and Bazley but with a higher value of air resistance. 

A much better agreement with data was obtained. The comparison is shown also 
in figure 4. 

The activities in ground-impedance measurement so far have generated suf- 
ficient data to establish trends for expected values of ground impedance. How- 
ever, the data are insufficient for rating all surfaces commonly encountered in 
aircraft noise measurement. For example, none of the ground conditions in the 
previously mentioned data group can adequately match the semidesert country of 
the western United States where many flyover tests are performed. Most ground- 
impedance data collected so far tend to follow the same trend as those summa- 
rized by Embleton et al. (ref. 2). In this type of data, both the real and 
imaginary components of the impedance are decreasing functions of frequency. 

In contrast, the data obtained by Dickinson show that the real component of the 
impedance remains approximately constant within the frequency range of 250 to 
1000 Hz. 
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One recognized difficulty is that existing methods of ground-impedance 
measurement produce data with a large amount of scatter. Part of the data 
scatter is a result of natural inhomogeneity of the ground -surface properties. 
Better methods for ground-impedance measurement are therefore needed. Theoret- 
ical modeling of ground impedance can also help in understanding the behavior 
of ground impedance. In addition to the studies by Delany and Bazley (ref. 28) 
and Chessell (ref. 39), the recent study by Donato (ref. ^0) contains theoreti- 
cal results comparable to the measured data. The development of theoretical 
models of ground impedance is important since it may lead to alternate methods 
of ground-impedance measurement. 


STANDARD PRACTICE IN CORRECTING FOR GROUND EFFECTS 

Most current standard methods for the correction and prediction of ground 
effects are empirical. In the area of long-distance propagation, simple 
methods of prediction are commonly used. (See ref. 1.) Sound attenuation is 
given as a function of elevation angle in these methods. However, some confu- 
sion exists regarding the range of elevation angle within which ground effects 
are important . According to recent theoretical results , the observed phenomena 
can be described analytically in a straightforward manner. The dependence of 
ground effects on elevation angle and ground impedance can be computed by using 
equations in the recommended prediction method and in the appendix. Much 
attention is given to methods of data correction for ground-based engine-noise 
and flyover-noise measurements. In this section, common practices in the lat- 
ter category are discussed. 

The objective of data correction is to recover the free-field spectrum 
from measurements with ground effects. The two possible approaches are either 
to minimize the acoustic problems associated with reflection or to correct for 
the ground effects during data processing. 


Modification of Physical Environment 

Flush-mounted microphones . - The first dip caused by multipath interference 
over a rigid surface occurs at a frequency given by 


r-|C 



where r-| is the source-to-microphone distance, c is the speed of sound, 
and h and z are the heights of source and receiver, respectively. By reduc- 
ing the microphone height, fo oan be made as high as necessary. For a dis- 
tance of 60 m, for example, and a sound source 2 m above the ground, fg will 
be above 10 kHz for microphone heights below 25 cm. Since the direct and the 
reflected waves will arrive at the microphone approximately in phase, the sound 
pressure is double the free-field condition. Therefore, 6 dB should be sub- 
tracted from all frequency bands for free-field equivalence. 
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Flush-mounted microphones are now being extensively used for outdoor 
acoustic measurements, especially for engine static tests. Owing to the 
effects of ground and wind and temperature gradients near the ground, precau- 
tions should be taken to guard against high-frequency losses in this configura- 
tion. Data above 1 kHz should be compared with equivalent bands as measured on 
a microphone at least 1 m above the surface. It is common practice to combine 
the low-frequency data from flush-mounted microphones with high-frequency data 
from the raised microphones, thereby reducing the flush-mounted-microphone data 
by 6 dB and the raised -microphone data by about 2 dB to obtain a composite 
free-field sound pressure level spectrum. 

The flush-mounted microphone should be used only over hard surfaces since 
the phase shift for reflection from a porous surface can introduce serious fre- 
quency distortion in the spectrum. For aircraft flyover measurements, good 
results have been obtained by facing a sheet of plywood approximately 1 .25 m 
by 1.25 m in dimension with aluminum and mounting the microphone in a hole 
drilled in the center. Furthermore, flush-mounted microphones should not be 
used for sideline noise measurements since ground effects, such as the surface 
wave, can introduce errors in the low-frequency data. 

Raised microphones .- It is practical also to arrange microphones high 
above the ground surface so that the first Interference dip falls below the 
frequency range of common interest. For example, a separation distance of 20 m 
combined with source and microphone heights at 6 m will bring the first inter- 
ference dip in the spectrum to about 50 Hz. Presumably the one -third -octave 
band at 50 Hz is not considered significant , and the higher order interference 
patterns are narrow compared with the one -third -octave bandwidth and are thus 
less influential on the measured data. This method is susceptible to errors 
caused by environmental factors and appears to be less accurate than the flush- 
mounted-microphone concept . 

Gravel and similar partly control led Impedance surfaces .- Many aerospace 
organizations use gravel or sand areas for static acoustic measurements to 
avoid the cost of large concrete pads. This practice gives reasonable repeat- 
ability of measurements for microphones which are at least a meter from the 
surface. The acoustic impedance of gravel test areas has been found to be 
relatively constant provided that they have sufficient thickness, are well 
drained, and remain above the water table. Since the engine static test facil- 
ities have a high rate of utilization, the gravel-filled area is sometimes 
acoustically calibrated for ground effects . The effects of the wind vector are 
often included in the calibration. 


Numerical Ground -Effects Correction 

Data correc tion with known ground impedance .- If reliable data of ground 
Impedance are known, ground effects in measured data can be corrected by means 
of analytical calculations. Figure 5 shows a comparison between a one-third- 
octave band spectrum of a jet engine on static test as recorded over a gravel 
surface with a microphone 75 m away and as recorded from a balloon the same 
distance above the source. The first microphone and the sound source were 4 m 
above the gravel surface. The presumed ground reflection characteristic is 
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determined by the difference between the two spectral curves and is shown in 
the lower half of figure 5 as a solid-line curve. The points are the calcu- 
lated values of the ground reflection characteristics using an impedance value 
taken from Dickinson's measurement for a broken tarmacadam surface. Even with 
the rough impedance estimate, there is a very good agreement with the measured 
data. Other discussions of data correction with known acoustic properties of 
the ground are given in references 41 and 42. 

Miscellaneous techniques .- Data which show prominent dips in lower fre- 
quency bands are frequently corrected by the trained eye and hand of the test 
engineer. This correction can be fairly successful for data obtained over hard 
surfaces if the test engineer understands the theoretical background and if the 
expected free- field spectrum is guided by, for example, model test data obtained 
in an anechoic chamber. In common practice the peak in the lowest frequencies 
will be lowered by 3 to 4 dB, a smooth line will be drawn through the first dip, 
and the high-frequency bands will be lowered 1 or 2 dB with peaks and troughs 
smoothed out. This procedure will obviously be incorrect in cases where a 
strong tone is present. Therefore, the test engineer should also be familiar 
with the characteristics of the sound source so as to distinguish which fre- 
quency bands legitimately should contain interference corrections and which 
ones may include real tones. He is also responsible for recognizing erratic 
behavior of the data, such as electronic noise sources. 

Computer algorithms performing the same sequence of operation have 
appeared in recent years. This procedure has been used by research organiza- 
tions where high-volume data processing is required. The computer-based method 
can utilize the theoretical results for ground effects in a more positive man- 
ner than the manual method. 

Note also that correction by visual smoothing can lead to serious error 
for data taken at long distances over a surface with finite impedance. The 
high-frequency acoustic signal can be systematically reduced through the action 
of ground effects. Consequently the ground-effects correction may require an 
increase from the measured level instead of the more common subtraction of 1 
or 2 dB. 

Analytic correction method s.- A number of computer-based methods for the 
correction of ground effects exist. Each method is designed for use at a spe- 
cific facility. However, these methods are definitely applicable at other 
facilities of similar type. The method given by Miles (refs. 43 and 44) is 
designed for an asphalt-surfaced acoustic test area where the typical distance 
of measurement is 20 m and the source and microphone heights are approximately 
4 m above the ground. This method is an iterative scheme where a model free- 
field spectrum is postulated at the beginning of the iteration. Ground effects 
are added to this spectrum and the result is compared with the measured data. 

A special error function is chosen to measure the agreement. If the error 
limit is exceeded, the parameters for the model spectrum and the ground effects 
are modified and a second iteration will be completed. Owing to the relatively 
short distance of measurement, a multiple-source location problem is often 
encountered for large test objects at this facility. This method has been 
found to provide satisfactory results for ground-effects correction under such 
conditions. 
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II 


Another approach is often adopted for ground -effects correction in gravel- 
filled test areas with typical source-to-microphone distances of 40 m to 60 m. 
The facility is calibrated acoustically to determine some typical parameters 
of ground impedance and expected ground effects on acoustic measurements. The 
computer program initiates a search for the frequency and magnitude of the 
first interference dip and then applies the appropriate ground-effects correc- 
tion for the entire sound pressure level spectrum according to the information 
obtained in the acoustic calibration of the facility. These methods are also 
found to be successful. 

Cepstral technique .- Another useful technique for ground-effects correc- 
tion was developed by Miles et al. (ref. 45). For acoustic measurements at 
short range with the sound source and the microphone high above a hard surface, 
the maximums and minimums of the interference pattern occur at regular frequency 
intervals. Furthermore, the first interference dip often occurs at a very low 
frequency which lies below the range of practical interest so that it may be 
ignored. A Fourier transformation of the logarithm of the spectrum results in 
a cepstral function in which interference maximums and minimums are represented 
as a single spike. The location of this spike is determined by the frequency 
span between two consecutive maximums or minimums. By removing this spike from 
the cepstral function and taking an inverse Fourier transformation, the result 
is a logarithm spectrum without ground effects. Note that this technique is 
applicable only for acoustic measurements over a hard surface where the surface 
reflection does not introduce any significant phase shift into the acoustic 
signal. 

The importance of ground effects has been recognized, and effective mea- 
sures have been taken in major test facilities for the inprovement of physical 
environments for outdoor acoustic measurements. In the area of data analysis, 
the methods now in existence show a trend of continued inprovement. In partic- 
ular, successful computer algorithms of ground -effects correction are now avail- 
able for specific conditions of test environment. The next logical step is to 
implement fully the advanced theoretical results in a computerized scheme. 


RECOMMENDED PREDICTION METHOD 

Specific formulas, obtained from recent theory, are suggested herein for 
the purpose of making predictions of the acoustic field in real experimental 
situations. These formulas are, as previously mentioned, based on the approxi- 
mate results of Chien and Soroka (ref. 11) and supplemented by the exact inte- 
gral solution of Thomasson (ref. 12). The physical model of sound propagation 
over ground which forms the basis of these theories has been discussed in the 
"Introduction.” The essential features of this model are the assumptions of a 
flat, locally reacting ground surface and a uniform atmosphere, and to the 
extent that these assunptions are justified for the particular experiment in 
question, the formulas can be expected to yield accurate predictions. Note 
that these formulas were derived for the case of time-harmonic waves radiated 
by a point source and hence can be applied directly only to this case. How- 
ever since these formulas represent Green's function solutions, they can be 
used to treat more general situations by means of superposition. For example, 
the solution for the case of a distributed source can be obtained by integrating 
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the point-source solution, weighted by the appropriate source term, over the 
source volume. A similar approach can be used to treat the case in which the 
source has a continuous frequency spectrum. 

The geometry of the situation is shown in figure 1 . A point time-harmonic 
source S is located at a height h above an impedance boundary which is 
taken to be coincident with the x,y plane. The distances from S and S' 

(the image source) to the receiver R are denoted by r-| and r 2 , respec- 
tively; r is the horizontal distance between S and R; and 0 is the angle 
of incidence. 

In mathematical terms, the problem of interest involves the solution of 
the equation 

(V2 + k2)<t> = 6(x)6(y)6(z - h) (1) 

in the region z > 0 of x,y,z space subject to the boundary condition 
3 (|) 

— + ikV(j5 = 0 (2) 

3z 


on z = 0. Here is the (complex) acoustic velocity potential and k = co/c, 
where w is the circular frequency (assumed positive) and c is the speed of 
sound. Also v is the specific boundary admittance, which is assumed con- 
stant, and is written in the form v = v-| + ±V 2 where Vi ^ 0. In addition 
to the boundary condition given by equation (2), a radiation condition corre- 
sponding to outgoing waves is understood to be imposed at infinity but is not 
explicitly indicated here. A time-harmonic factor common to both the 

source and the wave field, is also understood. The complex acoustic pressure 
amplitude p is given in terms of the velocity potential by the relation 
p = -ia)pQ((), where Pq is the unperturbed density of the acoustic medium. 


As noted previously, the various approximate solutions obtained by Chien 
and Soroka cover a wide range of practical situations, provided only that 
kr 2 » 1. Within this rather mild constraint, several subcases appear natu- 
rally and are defined below along with the corresponding approximate solutions 
(written in the present notation), which are recommended herein for prediction 
purposes . 


(1) ]v| « 1. In this case, the appropriate expression for the wave field 

is obtained from equation (25) of reference 11, and can be written in the form 


<|) = - 


1 


ikr-| 


ikr^ 


+ [r + (1 - D F(a)}- 


4tt 1 ri ~ ~ ^2 

where T is the plane wave reflection coefficient, defined by 


(3) 


cos 6 - V 
cos 0 + V 


( 4 ) 
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The function F is given by 


F(a) = 1 - erfc(o) 

where a is defined by 

/ 0/2 
/krsX 

a = j (oos 0 + v) 


( 5 ) 


( 6 ) 


In equation (6) the principal branch of the square root is understood; that is, 
the square root of any complex number ^ is defined by 

^1/2 = l^p/2 ®*P^~ i ^ 

-IT IT 

where -it < arg ? ^ it. Since — < arg v ^ , the argument of o lies in the 

2 2 

-3lT IT 

range i arg a $ — . The solution given by equations (3) to (6) has essen- 

4 4 

tially the same form as those of Ingard (ref. 8, eq. (13)) and Lawhead and 
Rudnick (ref. 9, eq. (21)). These authors did not, however, make clear the 
restrictions on the validity of their solutions. 

The power-series expansion 


20 


erfc(a) = 1 — = e~®^ ^ 

^ ^ (2n + 1)!! 


r2n+1 


(7) 


of the complementary error function vdien inserted into equation (5) yields the 
expansion 


F(o) = 1 - NlTrae*^^ + 2 o2 


(2o2)" 


^ (2n + 1)!! 


( 8 ) 


which is convenient for calculations when jol « 1. Recall that 
(2n + 1)!! = 1 • 3 • 5 • ... • (2n + 1) 
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Similarly, the asymptotic expansion 


.-a2r 


erfc(cr) s 2U(-Re a) + 




TTO 


1 - 


(2o2)^ 


(9) 


leads to the result 


F(a) s -2(/¥u(-Re a)ae'^‘^ 


2a2 (2a2)^ 


+ • 


( 10 ) 


which is useful when lcr| » 1. Here U denotes the unit step function, which 
is defined by 


r. 


u(s) = < - 

' 2 




(s > 0) 


(s = 0) 


(s < 0) 


The surface-wave term is easily verified to be implicit in the approximations 
given by equations (8) and (10) by simply truncating both expansions after two 
terms and inserting the result into equation (3). For the special case in 


which 6 = — (i.e., both source and receiver on the boundary), this procedure 

2 

yields 


jikr 


. 1/2 


4) = - 


+ - V 


(,.|v2)kr-l] 


2irr 


2 V27rr. 


( 11 ) 


when Id] « 1, and 


(j) = - 


le 


ikr 


2iTv2kr2 


+ U(-Re a)v 


.1/2 


l27rr; 


( 12 ) 


when ]a) » 1. The second terms on the right-hand sides of equations (11) 
and (12) can be identified as surface-wave terms. Note that the surface-wave 
term in equation (12) appears only when Re a ^ 0, or what is equivalent, when 

IT TT 

- — S arg V S — . 

2 4 
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Note that, although the original derivation of equation (25) of Chien and 
Soroka (ref. 11), or, equivalently, equation (3) of the present paper, required 
the additional assumptions that kr » 1 and k(h + z)^/r « 1, Chien and 
Soroka (ref. 11, sec. 6) showed that these additional conditions can be 
relaxed. 

(2) |v| > 1 (i.e., jv] is of order 1 or greater). In this case the 

appropriate expression for 4> is given by equation (12) of reference 11, 
which, after some manipulation, can be written 


(1> = $ 

4) = $ + 'i' 


(Vg ^ -f(Vi,6)) 
(V2 < -f(Vi,e)) 


(13a) 

(13b) 


where 


$ = — ■ 
4it 


feikr^ 

2iv 1 + V cos 9 

V 1 

eikr2) 

1-1 

i + 

^^2 (cos 9 + v)3 

r2 J 


(14) 


'1' = — kve” ivk(h+z)H^‘')[(1 _ v2)‘'/2ijr] 


(15) 


and 


2 - 1/2 

f(v-|,9) = CSC 6 (cos 6 + v-|)(1 + v-| cos S)(1 + 2v-| cos 6 + Vi ) 


( 16 ) 


Here '1' is the surface-wave term, and Hg denotes the Hankel function. The 
function f(vi,6) for selected values of 0 is plotted in figure 6. The 
surface-wave term in equation (13b) turns out to be negligible at all except 
near-grazing angles of incidence. To see this, note that, from figure 1, 
h + z = r 2 cos 9; also, as a consequence of the condition on v for the exis- 
tence of the surface wave and the assumption that kr 2 » 1 , note that -\>2 ~ 1 
whenever the surface wave exists. Hence 

|e-ivk(h+z)] < g-kr2Cos9 


which, since kr 2 » 1 , is negligible whenever cos 0 s 1 . Note also that the 
condition jv] ~ 1 of this subsection can be relaxed at all except near- 
grazing angles of incidence. This follows from the fact that the actual condi- 
tion for the validity of equation (12) of reference 11 is |cos 0 + v] 2 kr 2 >> 1. 
Therefore, provided that kr 2 » 1, the assumption |v| ~ 1 is needed only when 
cos 0 « 1 . 

Up to this point, the results of this section require for their validity 
that kr 2 » 1 . In order to treat the case in which kr 2 ~ 1 , it is necessary 
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to resort to numerical evaluation of the exact integral expression for the wave 
field. For this purpose, the integral expression of Thomasson (ref. 12, sec. V) 
appears to be most suitable. This solution can also be used as a check on the 
approximate results given previously in this section. 


In the present notation, Thomasson *s solution can be written 


ikri ikr2 


4it\ ri 


+ H 


r2 


where 


H = F 

H = F + G 


(V2 ^ -f(Vi,0)) 

(V2 < -f(Vi ,9)) 


(17) 


(18a) 

(18b) 


Here f(v-|,9) is given by equation (16) and 


F 


2lT 


(19) 


where 


and 


1= dt 

•^0 


D = (cos 0 + v)2 + 2i(1 + V cos 6)t - t^ 


( 20 ) 


( 21 ) 


In calculating the square root of D, the principal branch, as defined previ- 
ously, is to be used except when \>2 < -f(v-|,6) and t > t-| , where 

ti = -V 2 (cos 9 + Vi)(1 + v-| cos 0)“^ 


in which case the negative of the principal value is to be used. The quantity 
G in equation (I8b) is the surface-wave term; that is. 


G = - kVe-i^^(h+z)H^'>)[(1 . v2)^/^kr] 
2 


Note that the 
same as given 


surface-wave term and the condition for its existence are the 
by equations (13b) and (15). 
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PROBLEMS IN APPLICATION 

The preceding theoretical development is based on a set of idealized con- 
ditions which are not always met in practical applications. Some common prob- 
lems in this direction are now discussed. 


Finite Source 

The effects of finite-source dimension have been analyzed in detail by 
Thomas (ref. 46). In a fixed geometrical configuration where both the source 
and the microphone positions are high above the ground surface, ground effects 
on sound propagation are a result of phase differences between the direct and 
the reflected waves caused by ground reflection and path length differences. 

The maximum destructive interference occurs at a specific frequency where the 
total phase mismatch is 180^. However, such frequency is sharply defined only 
if the sound is emitted from a point source. In practice, a single jet engine 
may have a diameter greater than 2.5 m. 

Consider an engine noise test at static condition. The center line of the 
engine and the microphone are both positioned at 5.5 m above the ground with a 
horizontal separation distance of 60 m. Furthermore, assume that the ground 
surface is paved with concrete. In this configuration, the path difference 
computed from the center of the engine to the position of the microphone is 
approximately 1 m, and the corresponding frequency for maximum cancellation is 
approximately 174 Hz. If the source positions along the lip line at the engine 
intake are considered, the maximum cancellation frequency ranges from 141 Hz 
for points farthest from the ground to 224 Hz for points nearest to the ground. 
Hence, the acoustic energy emitted near the lip line will have a range of can- 
cellation frequencies extending over at least eight-tenths of an octave. The 
measured acoustic spectrum will show a broad but shallow dip in this frequency 
range, rather than a sharp minimtim. 


Filter Bandwidth 

The modification of the interference pattern by finite filter bandwidths 
has been analyzed by Howes (ref. 47) and discussed in reference 1. The general 
effect of the one-third-octave filters commonly used is to average out peaks 
and dips within bands above 1000 Hz and to broaden the lower ones similar to 
the effect of extended source dimension. The filters do not reduce much of the 
depth of the first interference dip since the width of this dip is larger than 
the filter bandwidth. 

Observe that the filtering procedure may change the measurements but has 
no effect on the physical phenomenon. Thus if there is a prominent tone at a 
single frequency in a high one-third-octave band, it may be increased or can- 
celed by multipath interference. The average ground -effects correction for the 
entire band would not be applicable for recovering the correct value of this 
tone. 
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Partial Coherence 


For noise emission with a broadband spectrum, such as jet noise, the 
origin of noise is often random in nature. The idea of phase cancellation 
should be treated with special care. The amplitude of the combined direct and 
reflected waves depends on the maintenance of coherent phase relations between 
wave components arriving at the point of measurement by different paths. The 
phase relation can also be lost through the effect of random perturbations 
along the path of propagation or as a result of reflection on a nonhomogeneous 
ground surface. The effect of coherence on sound pressure level can be shown 
in a simple analysis. A coherence coefficient can be defined as 

C(w) = Q(oj)/[So(w)St,(to)]^^^ (22) 


where Sq(o^) and Sj,(to) are the sound power spectral density of the direct 
and the reflected signals, respectively, and Q(w) is the magnitude of the 
cross-spectral density function between the direct and the reflected waves. 
(Note that this function is not always measurable in the configuration of 
ground reflection. It is well defined, however, when both the direct and 
reflected acoustic signals are known.) The coherence coefficient C(oj) can 
be considered as the fraction of initial acoustic energy in which phase rela- 
tion is maintained throughout the propagation process . The balance of the 
acoustic energy is distributed in waves of the same frequency with random phase 
relations . 

At the point of recombination of the direct- and re flee ted -wave compo- 
nents, the coherent portion of the wave will be summed according to the ampli- 
tude and phase of the two components , whereas the Incoherent portion of the 
wave energy will be summed according to the values of the acoustic intensity. 
This leads to the equation 

p2 

= C(03)[l + 2R cos (a + kd) + r 2] + [1 - C(o))](1 + r2) (23) 

PO^ 

where p represents the combined sound pressure, pQ is the expected free- 
field sound pressure at the microphone position, d is the length difference 
between the reflected path and the direct path, and R and a are the ampli- 
tude and phase of the reflection coefficient, that is, 

r(03) = R(to)ei“(“) (24) 

where r(w) is defined in equation (4). The values of R and a are fre- 
quency dependent. Since the coherence coefficient is not always a measurable 
function, some assumptions may be necessary for practical applications. A 
reasonable choice is a Gaussian distribution 

C(w) □ (25) 

where a is a real and positive constant. 
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At the frequency of maximum destructuve interference, the direct- and 
re fleeted -wave components are 180° out of phase at recombination. Equa- 
tion (23) can be simplified to give the relative sound pressure level in ref- 
erence to the expected free-field sound pressure level (SPL) at the microphone 
position 


ASPL = 10 logio{c(w)(1 - R)2 + [1 - C(w)](1 + R^)} 


(26) 


Calculated values of this interference dip for different values of ,C(w) and 
R are given in figure 7. Note that the presence of a small fraction of inco- 
herent wave can make drastic changes in the value of the interference attenua- 

tion. For example, with a reflection coefficient of 0.90, the magnitude of the 
interference dip reduces from 20 dB for a wave with perfect coherence to merely 
7.9 dB for a wave with a coherence coefficient of 0.95. 

In the case of jet noise, sound emitted in different directions will not 
be perfectly coherent. According to experimental indications, the coherence 
factor (ref. 48, figs. 20 and 21) is 

C(to) = cos'* 3 (27) 

where 3 is the angle between the direct and the reflected ray paths at the 

source. For an emission angle difference of 10°, which is common for typical 

engine static test configurations, the coherence factor is approximately 0.94. 
Furthermore, jet noise is coherent only over a finite length of time. The 
typical coherence time scale for a large jet engine is of the order of 10 to 
15 msec. For a propagation path difference of more than a few meters, the time 
incoherence of the acoustic signal must be taken into consideration. Such con- 
ditions may be encountered for flyover measurements with microphones positioned 
high above the ground. 

Under conditions of long-range propagation, several additional mechanisms 
may become significant in loss of signal coherence. Among the most important 
of these are the effects of atmospheric turbulence and the effects of random 
ground roughness. An analysis of the former has been given by Ingard and 
Maling (ref. 49); however, that analysis considered only single scattering and 
hence is valid only for relatively limited propagation distances. With regard 
to the latter, mathematical techniques which include multiple-scattering effects, 
and which are therefore applicable to long-range propagation, have only recently 
been developed. (See, e.g., refs. 50 to 52.) No attempt has yet been made to 
apply these techniques to real problems involving propagation of sound over 
ground. It thus appears that predictions of coherence loss due to these mech- 
anisms must remain, at least for now, largely empirical. 

An additional effect of atmospheric turbulence, which is not related to 
the presence of a nearby reflecting surface, is manifested in the appearance 
of random fluctuations in the sound pressure level. Observations indicate that 
these fluctuations increase initially with propagation distance, but that at 
longer distances they approach a limiting, or saturation, value. (See, e.g., 
ref. 3, sec. IIIC.) This saturation phenomenon, which is a general feature of 
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wave propagation in a random medium and vdiich is not amenable to a single- 
scattering analysis, has been the subject of a number of recent theoretical’ 
investigations. (See, e.g., ref. 53 and references cited therein.) 


Wind Gradients 

In the presence of wind, the analysis of ground effects becomes very dif- 
ficult from a theoretical point of view. However, in geometrical configura- 
tions typical of ground-based engine noise measurements, the effects of a small 
wind vector can be included by means of simplified analysis. The phase differ- 
ence frequency is determined by the phase shift caused by the ground reflection 
and the propagation path difference. In the presence of a wind vector, the 
difference in average propagation velocity between the direct and the reflected 
paths will produce an additional segment of path difference. For a typical 
boundary-layer profile, this additional path difference can be given as 

d-| H 0.3rv/c (28) 

where r is the horizontal distance, c is the speed of sound, and v is the 
wind-velocity component along the direction from the source to the receiver. 

In equation (28), the source and the microphone are assumed to be at the same 
height above the ground. According to recent studies, the deviation of the 
frequency for maximum interference attenuation can be approximately two-thirds 
octave on either side of the expected peak attenuation frequency. For practi- 
cal applications the effect of wind on the wave-propagation path difference 
should be included in the computations of ground effects. 


CONCLUDING REMARKS 

A unified method has been recommended for the prediction of ground effects. 
This method is applicable to long-range sound propagation, correction of data 
in engine test-stand configurations, and aircraft flyover noise measurements. 

The results given in this paper contain sufficient detail for the purpose of 
numerical calculation of ground effects. 

Sufficient data are available at this time to form a band of expected val- 
ues for ground impedance for common surfaces such as grassland, sand, and soil. 
Measurements from different sources seem to agree in general. However, there 
is a large variance of data scattering. In order to obtain more accurate data 
on ground impedance, methods of measurement should be improved. Furthermore, 
the measurement of ground impedance using existing methods is time consuming. 
Improved procedures may encourage researchers to make on-site measurements of 
ground impedance in conjunction with acoustic measurements. 

Theoretical analysis for applications to aircraft noise prediction is 
adequate for the present. However, the analysis of coupled effects of ground 
and atmospheric perturbations on sound propagation near the ground seems to 
provide the possibility of a better description of the natural environment. 
Further research in the area of layered representation of the ground surface 
may also be beneficial to practical applications. 
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From an analytical point of view and also from evidence shown in existing 
experimental results, the recommended method appears to be sufficiently accu- 
rate for practical applications. However, the method should be validated 
before it is accepted for general industrial applications. Since improved 
instrumentation and facilities of data analysis are now available, an experi- 
ment for sound propagation over the ground with careful documentation of ground 
impedance and atmospheric conditions should be conducted. 

A validation program is beneficial from another point of view. In some 
field measurements, the ground-impedance environment will not be known beyond 
a general description of the physical appearance of the ground condition. In 
the process of validation, experience can be gained so that a band of nominal 
values for various parameters can be recommended. Thus, an estimate of ground 
effects can be obtained in the absence of direct information concerning the 
environmental factor. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 8 , 1 977 


25 



APPENDIX 


RELATIVE IMPORTANCE OF PLANE WAVE APPROXIMATION 


A general expression for the sound pressure field near an impedance bound- 
ary can be given as 


P = 


eikr^ 

4uri 


{r(6) + [1 - r(6)]F(a)} 


4iTr2 


(A1) 


Symbols in equation (A1) are defined the same as those in the main text. If 
the effect of the surface wave is not considered, then 


[1 - r(6)]F(a) = 


2iv 

kr 2 (G.os 0 + v)3 


(A2) 


when ] kr2\)2| » 1 and v cos 0 « 1 . For sound propagation over a long hori- 
zontal distance, the path-length difference between the direct and the reflected 
waves is very small. Therefore 

r z r-i z P 2 (A3) 

and the direct and the reflected waves can be combined into one term (ref. 31) 

gikr 

P 1 = cos 0 (A4) 

2TTr(cos 0 + v) 


IT 

An angle 0 * can be defined such that 0 ' = — - 0 . For a receiver height that 

2 

is small compared to the source height, or conversely, the angle 0' is 
approximately equal to the elevation angle of the source. For small values of 
the angle 0 ' 


kh 

cos 0 = sin 0* =0* s — (A5) 

kr 

By substituting equations (A2) to (A5) into equation (A1), the sound pressure 
field can be given in a simple form 

eikr 

p = ^ 

2iTkr'^(cos 6 + v)v 

It is clear from equation (A6) that the combination of direct and reflected 
waves given by equation (A4) dominates when 


khv + it 


iCOS 0 + V 


(A6) 
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khvl » 1 


APPENDIX 


(A7) 


This result leads to a number of Interesting conclusions as follows: 

(1) Both the combination term and the ground-wave term (according to the 
definition of ref. 27) attenuate at the rate of r"^. These two terms are not 
separable in measurements at long distances away from the sound source. 

(2) If the condition specified in equation (A7) is met, the ground-wave 
term can be neglected. The overall sound pressure is directly proportional to 
the elevation angle only when cos 6 « Iv] . (See eq. (A4).) The ground 
effect is negligible for cos 6 » jv] . Therefore, this equation provides an 
analytical basis for the empirical methods for estimating ground attenuation 

by using the elevation angle alone. For aircraft in-flight sideline noise mea- 
surements, kh is normally large and equation (A7) is satisfied. However, the 
ground attenuation should be considered as a function of frequency since the 
ground admittance is a function of frequency in most cases. 

(3) This analysis also provides a limit to the validity of the previously 
mentioned empirical method. It should not be applied to the low-frequency 
range or under conditions vrtiere the surface wave may have a strong influence. 

Finally, note that Piercy et al. (ref. 31) pointed out that the ground 
wave is a correction of the sphericity of the wave front. Naturally, the 
ground-wave term becomes less important at distances far from the source since 
the sound field can be better approximated by a plane wave. 
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Figure 2.- Measured acoustic impedance of sand and grassland. 
(From Dickinson and Doak, ref. 32.) 
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Figure 3.- Magnitude of acoustic impedance of grassland, 
chipped granite, and gravel. 
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Figure 4.- Comparison of estimated acoustic impedance of ground by Delany 
and Bazley (refs. 28 and 29) and Chessell (ref. 39) with measured data 
reported by Embleton et al. (ref. 2). 
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